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Contract  A?  29(601) -4136 
Seal-Annual  Technical  Summary 
for  the  period  ending  31  One.  1962 


The  work  discussed  In  the  laat  seal-annual  sunmary  Is  pro¬ 
gressing  toward  completion  and  results  are  being  prepared  for  Inclusion 
In  the  final  report. 

In  addition  to  the  work  described  there,  a  machine  program  for 
solution  of  the  spherically  symmetric  coupled  radiative  transport  and 
hydrodynamic  equations  has  been  formulated. 
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Radiative,  Cooling  at  Very  Late  Times 

To  provide  a  means  of  calculating  radiative  cooling  at  very 
late  times,  such  that  hydrodynamic  motions  are  relatively  unimportant^^ 
and  the  hot  gas  is  everywhere  optically  thin,  a  simplified  code  is  being 
developed  to  calculate  the  rate  of  radiative  cooling  for  gases  at  temper¬ 
atures  below  10,000  subjected  to  a  given  variation  of  pressure  versus  time. 

•  *  c 
For  this  calculation  the  ealssivitles  of  Kivel  and  BaileyX  are  employed , 

together  with  newly  developed  analytic  fits  to  the  Gilmore  tables^  This 

code  is  expected  to  be  most  useful  for  low  altitude  detonations,  where 

the  late-time  motion  can  be  reasonably  approximated  by  simple  models. 
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'  Central  Description 

The  code  developed  for  late  time  fireball  computations  la  a  gen¬ 
eralisation  of  the  Ricbtmyer-von  Neumann  artificial  viscosity  technique  for 
a  spherical  Lagranglan  system  in  which  the  divergence  of  the  radiative  flux 
appears  in  the  hyperbolic  partial  differential  equation  which  expresses  the 

conservation  of  energy  for  a  given  mass  sons*  The  basic  differential 

(2) 

equations  are  those  given  by  Rlchtmyer, 
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modified  to  include  this  additional  radiative  term.  Without  this  term, 
the  solution  of  the  equations  is  a  standard  initial  value  problem  which 
can  be  integrated  explicitly.  Such  solutions  are  well  known;  only  those 
modifications  required  to  include  the  radiative  transport  will  be  discussed 
here. 

Formulation  of  the  Radiative  Energy  Transport 

To  understand  the  basic  problem  of  radiative  transport  consider  the 
standard  equation  of  radiative -transfer  for  a  media  in  local  thermodynamic 
equilibrium  and  for  a  plana  parallel  geometry: 

toiQ  =  KvfaCn*)  -  !,(*<**>) 

where  the  quantities  are  defined  in  the  standard  terminology  as  given  by 
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This  equation  represents  the  change  la  the  specific  Intensity  of 
the  radiation  along  a  ray  having  a  given  direction  in  space.  Since  an 
infinite  number  of  directions  for  rays  can  be  assigned  to  any  unit  volume 
under  consideration,  it  is  clear  the  eq.  h  represente  an  infinite  set  of 
differential  equations  which  are  functions of  position,  direction  in  space,  ' 
and  radiation  frequency. 

The  solution  for  this  set  of  equations  has  been 'obtained  in 'dosed 
'  •  (ibid  pages  55  ff)  1 

form  by  Chandrasekhar  based  on  a  method  in  which  the 

radiation  field  is  represented  by  2X  streams  where,  in  the  limit,  X 

,1 

approaches  infinity.  The  convergence  of  this  solution  is  shown  to  be 
critically  dependent  on  the  proper  choice  of  radiation  streams  usea  in 
the  first  approximation,  (l.e.  N  ■  l),  which  corresponds  to  the  two  stream 
method  of  A.  Schuster  and  X.  Scbwarxscbild  (1905).  Suitable  choice  of  the 
directions  of  the  tvp  rays  used  to  represent  the  radiation  field  results 
in  a  formulation  of  the  gross  energy  transport  whose  accuracy  is  considerably 
greater  than  the  accuracy  with  which  we  know  the  appropriate  absorption  co¬ 
efficients.  The  accuracy  of  this  approach  improves  when  it  is  applied  to 
thin  isothermal  slabs  as  contrasted  with  the  usual  application  to  an  atmosphere 
of  essentially  infinite  extent. 

The  basic  physics  of  the  Schuster-Schwarsschild  method  has  been  carried 
over  to  spherical  geometry.  Xn  this  application,  the  inclination  of  the 
chosen  radiation  streams  changes  with  the  relative  radius  of  curvature  of 
the  spherical  tones.  When  the  radius  of  curvature  becomes  large  compared 
to  the  sons  thickness  the  appropriate  inclination  is  that  which  is  applicable 
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for  plane  parallel  geometry.  Indeed,  la  tbit  Halt,  the  exact  solution 
Tor  the  flux  in  terns  of  exponential  integrals  any  be  employed  and  no 
approximation  is  required.  Kovover,  in  tbe  present  application,  other 
conditions  of  the  problem  made  it  desirable  to  maintain  tbe  two  stream 
approximation,  tbe  choice  of  inclination  being  determined  by  tbe  condition 
that  the  difference  la  these  tvo  solutions  be  held  to  a  minimum.  Tbe  proper 
choice  of  inclination  near  the  center  of  tbe  sphere  has  been  temporarily 
postponed  by  choosing  tbe  central  zones  such  that  they  are  optically  thick 
and  by  requiring  that  tbe  tvo  stream  solution  transform  into  tbe  diffusion 
approximation  which  is  applicable  under  these  conditions. 

Tbe  actual  formulation  of  tbe  flux  is  then  carried  out  using  tbe  formal 

<3); 

integral  solution  '  of  tbe  equation  of  transfer  for  tbe  specific  intensity 
along  tbe  representative  rays.  Tbe  source  function  is  assumed  to  be  that 
which  is  appropriate  under  conditions  of  local  thermodynamic  equilibrium 
and  is  expanded  in  terms  of  the  temperatures,’  and  derivatives  of  tbe  tem¬ 
peratures,  at  tbe  zone  boundaries.  While  tbe  code  has  been  written  to  allow 
tbe  spectrum  to  be  divided  into  any  desired  number  of  frequency  bands,  tbe 
initial  computations  have  been  carried  out  using  tbe  Bosseland  mean  absorption 
coefficients  for  tbe  purposes  to  be  enumerated  later. 

The  numerical  solution  is  carried  out  by  explicit  integration  of  tbe 
basic  equations,  except  for  the  energy  equation,  using  tbe  well  known 
Rlcbtayer-von  Neumann  method  Tbe  integration  of  tbe  energy  equation 
cannot  be  carried  out  explicitly  since  tbe  radiative  flux  term  for  any 
particular  zone  is  dependent  upon  tbe  temperatures  of  tbe  neighboring  zones. 
Since  tbe  density  is  known  at  this  stage  of  tbe  solution,  it  is  convenient 


y 


to  express  tbo  energy  equation  in  tens  of  none  temperatures  alone.  The 
solution  is  then  carried  out  by  assuming  an  approximate  temperature  for 
each  tone  and  using  a  Newton-Raphson  Iteration  scheme  to  adjust  the  sons 
temperatures  until  changes  in  temperatures  on  subsequent  iterations  are 
less  than  a  pre -determined  amount. 

This  method  of  solution  involves  the  solution  off  X  linear  equations 

a 

in  X  unknowns,  where  X  is  tbs  number  of  sones  which  is  usually  of  the  order 
of  a  hundred.  Such  a  system  is  in  principle  easily  soluble  by  computer 
methods,  but  the  procedure-:  is  time  consuming.  In  practice  only  a  fev 
zones  in  each  direction  measurably  Influence  the  temperature  of  a  particular 
zone  which  results  in  a  simplification.  In  a  first  solution  of  the  problem 
a  five  zone  centered  system  has  been  used  which  results  in  a  matrix  having 
five  terms  across  the  diagonal.  After  suitable  scaling,  such  a  system  is 
rapidly  soluble  by  direct  elimination  and  back  substitution. 

As  part  of  the  check  out  procedure,  the  equation  of  state  and  formulation 

(ibid) 

of  the  mean  absorption  coefficient  given  by  Brode  were  used  so  as  to 

reproduce  his  results  over  an  appropriate  period  of  fireball  history.  Using 

* 

100  sones,  and  3  iterations  per  hydrodynamics!  cycle,  the  code  proceeds 
successfully  about  200  time  steps  per  hour.  With  minor  modifications,  the 
code  has  also  been  applied  to  a  problem  in  thermal  ablation  in  fireballs  and 
to  problems  in  stellar  stability. 

Interpolation  formulae  are  currently  being  written  for  mean  absorption 
coefficients  for  chosen  spectral  bands  appropriate  to  late  fireball  develop¬ 
ment. 
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Initial  Energy  Deposition  in  Variable  Density  Atmosphere 


When  a  nuclear  device  Is  detonated  la  the  atmosphere ,  a  large 
percentage  of  the  energy  released  is  radiated  away  and  is  deposited  la 
the  surrounding  atmosphere .  Since  the  air  density  exhibits  a  large 
variation  with  altitude,  It  follows  that  the  scale  of  distances  over 
which  this  energy  is  absorbed  by  the  air  will  be  dependent  upon  the  burst 
altitude.  Consequently  the  source  density  say  also  exhibit  a  dependence 
on  altitude. 

The  sequence  of  radiative  and  hydrodynamic  events  which  follow  this  _ 
Initial  phase  are  dependent  upon  the  radial  energy  distribution.  Thus, 
one  would  expect  the  scaling  laws,  for  various  weapons  effects  as  a 
function  of  altitude,  to  also  exhibit  a  dependence  on  the  Initial  energy 
deposit.  Zt  Is  therefore  of  Interest  to  perform  computations  of  this  die* 

a 

trlbution  as  a  function  of  weapon  yield  and  altitude. 

Because  of  the  large  nuoer  of  cases  of  Interest,  a  machine  code  was 
written  to  perform  the  appropriate  computations.  The  bsislc  equations  In*. 
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where  ~  -  energy  per  unit  volume  deposited  at  a  distance  B  from  the  ’ 
point  of  detonation. 
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p  ■  air  density  as  a  function  of  altitude  A 
absorption  coefficient  of  air 

W y"  total  energy  of  frequency  '))  radiated  by  tbe  source.  Equation  A 
was  evaluated  using  an  empirics!  fit  to  tbe  ARDC  Standard  Atmosphere  (5). 

The  absorption  coefficients  used  were  those  of  Gilmore  (#}  for  cold 
air.  Stripping  vas  neglected  since  it  appears  to  have  only  a  snail  effect 
upon  tbe  later  development. 

In  coding  the  problem,  tbe  radiative  output  of  tbe  source  woe  divided 
into  11  separate  spectral  regions  and  a  corresponding  average  opacity  was 
chosen  for  each  of  these  regions.  A  minor  modification  of  the  code  allows 
a  larger  number  of  spectral  regions  to  be  used  when  desired.  Tbe  code  carries 
out  radial  computations  for  tbe  downward  direction  and  at  directions  in  space 
differing  by  10  degrees,  or  any  multiple  of  10  degrees,  as  may  be  desired. 

The  printed  output  for  each  direction  in  space  lists  the  radius  from  tbe 
point  of  detonation  for  35  different  points,  and  tbe  sir  density,  altitude, 
and  energy  density  for  each  point. 

Temperatures  were  determined  using  the  tabulated  equation  of  .state 

<7>  •  (8) 

data  given*  by  Gilmore  ,  end  Eilsenratb  end  Beckett  .  Initial  rune 
have  Shown  that  a  better  representation  of  tbe  equation  of  state  for  air 
is  desirable. 

Tbe  initial  runs  included  a  variety  of  yields  and  altitudes  including 
cases  applicable  to  previous  field  test  situations  and  cases  of  interest  for 
future  planning.  Considerable  success  was  obtained  In  predicting  a  number  • 
of  the  effects  observed  In  earlier  field  teste,  but  s  better  handl <"g  of  tbe 
equation  of  state  is  indicated  before  tbe  results  art  presented  in 


lav  form 
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Various  aspects  of  the  role  of  debris  is  determining  ionization 
levels  are  being  studied.  The  work  on  ionization  of  debris  atoms  has 
been  mentioned  in  previous  reports.  A  second  interesting  question  is 
bow  to  reconcile  the  large  dimensions  of  observed  debris  clouds  vitb 
the  small  values  found  in  hydrodynamic  calculations.  Finally,  we  note 
that  in  low  altitude  detonations,  dust-like  impurity  particles  may  be 
formed  and  play  a  significant  role  in  electron  recombination,  as  recently 
suggested  by  Rosen10. 
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